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1. INTRODUCTION 

At current days distributed generations systems (DGs) have been used to generate electricity 
because of their high reliability and low gas emissions, especially for CO2 gas, which made them 
environmentally friendly and low in pollution, and given the most important factor, which is cost, compared 
to the traditional source of electricity, it is more practical and less expensive, connected to DGs is a low 
voltage distribution tool for renewable power components which includes; Fuel cell, wind or photovoltaic, 
for the grid the usage of a power inverter [1]. DG units generate electricity closer to load centers, avoiding 
the expense of energy transportation and lowering transmission line power losses. Furthermore, when 
compared to centralized generation stations, DG solutions save more money [2]. As it is known that 
renewable energy sources generate direct current (DC) power, so it is necessary to use a power inverter to 
convert it into alternating current (AC) power, this is why it is used in these systems including active and 
reactive energy which can be controlled, However, in voltage source inverter (VSI), the pulse width 
modulation (PWM) approach causes the appearance of high and low harmonics, which are detrimental to the 
system's balance and must be eliminated [3]-[5]. That is why it's vital to utilize a filter to do away with these 
problems and treat them. In previous research, the inductor (L), inductor-capacitor (LC) and inductor- 
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capacitor-inductor (LCL)-filters are designed and they have the following characteristics. Despite its 
popularity and ease of use, a single inductor L-filter has a low attenuation and an excessive inductance rate. 
The voltage drop across the inductor causes negative machine dynamics, resulting in a long response time. 
The inverter switching frequency must have an excessive value to sufficiently attenuate the harmonics when 
using an L-filter [6]. L-type filters can attenuate -20 dB/dc, that's a small amount, so it needs a capacitor that 
produces a large and excessive impedance in the frequency variety. For modulations, the LC filter is ideal 
that have notably high load impedance and above running frequency. After that, it needs a large capacitance 
to reduce losses and thus cost. Nevertheless, this will also result in an increase in reactive power, In addition 
to the possibility of the appearance of the phenomenon of resonance [7]. Because it has a larger decoupling 
between the converter and the grid impedance, better attenuation of the switching frequency harmonicas well 
as reducing total harmonics distortion (THD), to compensate for the high capacitance, the LCL-filter can be 
utilized. Since this system's stability will be affected by the filter's resonance, passive damping is utilized by 
connecting resistors in parallel with the filter output capacitors to damper the resonant harmonics [8]. 

This work proposed the power control by using PI controller, which has a fast response to changes 
in the control process since it takes instant corrective action when an error occurs, and it has a very stable 
control process if the proportional gain (Kp) and integral gain (Ki) factors (parameters of PI controller) are 
met and properly selected. The PI controller will be using two methods: the first is the traditional trial and 
error approach and the particle swarm optimization (PSO) algorithm. Power injection controlled with sudden 
or gradual load changes, through the usage of an LCL filter between the inverter and the grid, the harmonic is 
removed and the THD is reduced to less than 5%. The technique of decoupling control is employed. The 
important characteristic should be iq=0, which gives the idea an easy configuration and secure manipu- 
late [9]. 

This search was completed in the following sections: in section 2, the model of LCL-filter obtained 
between converter and network. The equations for the system is obtained by using the direct decoupling 
method. In section 3 show the transfer function of system with LCL filter, section 4 represent stability 
analysis of LCL filter, section 5 power control strategies and section 6 represents the simulation results for 
each controller PI, finally, the conclusions drawn from this work are summarized in section 7. 


2. MATHEMATICAL MODEL OF LCL FILTER INVERTER 
The network's stationary frame analysis can be expressed as shown in Figure | using Kirchhoff's 
voltage law, as shown in: 


; d. 
vn-ni- ligu Ve =0 (1) 
d. ; 
Vv -y = Lgi Fry (2) 
z d. 
V = r2i2 — L2 ;zi2 72€ = 0 (3) 
d. ; 
V -e = Lz g t2 t "2i2 (4) 


where: iy = ic + iz, ig = CLV, , L = Lp + Lg , T2 = T + rg 

V,: The voltage of converter side. 

i, : The current of converter side. Ls: The inductance of grid side filter. 
V; : The voltage of capacitance. Lg: The inductance of grid. 

c: The current of capacitance. r: The resistance of converter side filter. 
i2: The current of grid. rg: The resistance of grid side filter. 

L, : The inductance of converter side filter. rz: The resistance of grid. 
The mathematical model inside the aß frame is: 


~. 


, d. 
Viag — Vcag = riliag + Li 1B (5) 
! d. 
Veag — Eag = T2iz2ag + Le 7 izap (6) 
, a 
leap = C zg Vcag (7) 
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itag = icag + i2aß (8) 
it is important to note that the two phase d-q system has a property, which is that any vector has an amount 


that rotates with the frame and with the same velocity, which is the angular velocity œ and the equations are 
as shown in: 


r d. > 
Via — Vea = Tilia + Ly ha — Liwia (9) 
; d. š 
Vig — Veq = filiq + Ly ha + L, Wiig (10) 
; d. ; 
Vea — ea = niza + Lz laa T LoWiza (11) 
s d. ; 
Voq = eq = T2l2q + Lz t24 + Lzwizq (12) 


P d r d 
where: led = CVa š leq = CVa 


The decoupling components are defined: i,g4°°°UP!e? = WLyigi > ig =o L, igi. So the currents 
will be: 
la = lca + tag + Wlyigy (13) 
lig = leq + lag — @Lr hia (14) 


The three phase voltages and supply currents at the point of common coupling (PCC) are converted 
into their respective components using Clarke's transformation. The stationary components are later changed 
into their equivalent rotational d and q components, synchronized to the supply voltage, using Park's 
transformation [10]. The phase locked loop (PLL) will trace the phase angle difference between the inverter 
and grid voltages and set it to zero, as well as a current controller to control the inverter's signals and measure 
the voltage magnitude that must be synced [11]. The use of a stationary or rotating frame is a typical method 
of analyzing three-phase systems Figure 2. 


DC/AC Inverter LCL-Filter Grid 


Figure 1. LCL-filter inverter connected to the grid 


Figure 2. shows a stationary and a rotary dq frame 
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3. TRANSFER FUNCTION OF SYSTEM WITH LCL FILTER 
Rearranging in (9) and (10) and putting S in place of: 


d ; A ; : 
at Viaq kan cdq + Yilidq + SLiliaq + Liwiq = Liwia 
ts i ; ; : 
Viaq = 5c leaq + Tiyaq + Slytag + Liw(isq T ia) (15) 
: a FAs 
where: Vegg = 5c beaa 
laq = leaq + loaq (16) 
In the same process, the carried out (11) and (12). 
Veaq — €aq = Tzizaq + SLzizaq + Low(ing — ina) (17) 
eaq = Veaq — Tai2aq — SLzizaq — L2W(i2g — ina) 
a : ; ; i 
Caq = zgicdq 7 "2İ2aq — SL2izaq — Lz%(izq — iza) (18) 


If the filter is perfect, harmonics higher than the fundamental harmonic do not appear in the LCL filter's 
output current (€gg = 0) [8]. 


icdq = SCTzingg + S?°CLzizaq + SCL2w(izq — ina) (19) 
Substituting as shown in (16) and (19) in (15) 


Viag = 5c |SCT2i2aq + S?CLoizag + sCLzw(izq — iza )| + Tilceaq + Mleag + SLiicaq 
+sLiizaq + L10(iq — bra). 


Viaq = Tzizaq + SLzizaq + L2%(izq — iza) + SCr1ů2izaq + S?CLatyizag + SCL2T (ing — ina) + 
Tiizaq + S?°LaCnzizaq + S°CLiLzizaq + S?°CLiL20(izq — ina) + SLiizaq + L10(izq — iia). 


Viaq = İzaq [S3CLy — 3? (Clar, + Clir, + Cla Lzw(izq — iza)) + 
s (Ly + La + Cryty + Cla @(ing — iza) ) + (r1 + r2 + L1o(iq — a) )] (20) 


ing x 
4 = —_____ , where: dg = +n , 
Vidq a3S°+a2S*+a,S+ag 


For (w = 0) So the transfer function become: 


a, = Li + Lz + CL,L, a2 = CLr + CLir , ag = CL,L, 


4. STEABILITY ANALYSIS OF LCL FILTER 

To minimize resonance harmonics, the passive damping method has been used, which involves 
adding passive elements to the filter circuit, resistors in series or parallel with the capacitor are commonly 
used [12]-[16]. The serial branches minimize losses while lowering high-order harmonic dampening 
capacity. The parallel department, on the other hand, has no effect on low and excessive-order harmonics, but 
it does increase losses. The adding resistor, on the other hand, will result in higher loss and less attenuation of 
high-frequency harmonics. Alternatively, active damping strategies are used to suppress resonance and 
stabilize the system, with the basic idea being to introduce more variables as damping terms to the present 
control loop [17]-[20]. The frequency responses are investigated in Figure 3 by substituting the following 
values of inductances, capacitors, and resistances in transfer function. Where: L4 = 5.3052 x 1074 H , L, = 
1.5321 x 1074 H , C = 1.10022 x 1074 F. 

Figure 3(a) shows the frequency response of closed-loop system when neglecting parasitic resistance 
of the filter inductance (i.e. %,72 = 0), while Figure 3(b) shows the frequency response of the actual 
inductance resistance rı = 0.1667 0,72 = 0.0096 Q. It’s noted from Figure 3 that when a certain value of 
resistors is added to the filter, the system becomes more stable than the purely inductive filter. 
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Bode Diagram 
Gm = -184 dB (at 8.74e+03 rad/s) , Pm = 90 deg (at 1.51e+03 rad/s) 
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Figure 3. The frequency response of the LCL-filter for, (a) r1, r2 = 0 and (b) % = 0.1667 0,7 = 0.0096 Q 


5. POWER CONTROL STRATEGIES 

Controlling the active and reactive components of the source through implementing control 
strategies on the inverter could be used to achieve power and voltage regulation techniques [21]-[25]. Grid 
converter power control is based on instantaneous power theory and, as a result, on a reference frame 
definition of power. The VOC is based on the rotation of a dq frame at a high speed with the d axis aligned 
with the grid voltage vector. The current id* is controlled to injected active power and iq* is controlled the 
reactive power. For making the grid current vector to be in phase with the grid voltage vector, iq* should be 
zero [8]. As a result, the grid converter's active and reactive power are shown in: 


ee : 

P= zeala + eqlq (21) 
3 ; 

Q= 5 Pata — Calg (22) 


where: 
P: Real active power, Q: Real reactive power 

Active and reactive power will be proportional to id and iq, respectively, supposing the d axis is 
correctly aligned with the grid voltage e, = 0. The block diagram of the power control system is shown in 


Figure 4. 
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P= egia (23) 


Q = Žeaiq (24) 


LoadL 


cB Load3 


Inverter | ——> LCL-filter 


$ PLL 
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Generator P 


Jf 
d dq 
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Figure 4. The control system's structure diagram 


5.1. PQ open-loop control 

Powers feed-forward control is the most straightforward ways to implement voltage-oriented 
control. The active and reactive power command signals are then translated into the reference current's d and 
q components using the matrix in: 


Alle as 


i [= Vga" +Vgq" ae 


where: 

Vja: Grid d-axis voltage, Vjqg : Grid q-axis voltage 

iq : d- axis required current, ig": q- axis required current, P*: Injected active power, 
Q*: Injected reactive power 


5.2. PQ closed-loop control 
The dynamics of powers control can be determined as a result of changes in grid voltage using 
closed-loop control: 


ia* = P*—P (26) 
ig =Q*-Q (27) 


where: P = Vaiga + Vgqigqa: Q = 


—Vjqia + Voat 
6. SIMULATION AND RESULTS 
6.1. PI controller with trial-and-error method 

To optimize the PI controller gain, many common tuning methods have been utilized, such as trial 
and error tuning methods to obtain optimal values for the PI controller values Kp and Ki. The PI controller 
was used to control the injection of the required power into this network, which is a generator that generates 
electrical energy and feeds three different loads that enter the work at certain times. This means that the 
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network generator will need to prepare a higher capacity in the event of high loads entering the work. The 
power from the renewable energy source will be injected into it to help it equip the required energy with the 
overall. System simulation shown in Figure 5. Table 1 shows how PI controller gains are chosen by trial and 
error. 
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— oe 
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22k/380V 250kVA 


Grid Generator 


‘Theos Phase Broa pater ig E 
Load 2 three-phase 
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Load 1 three-phase = 


balance load 


Figure 5. Simulation model of power control with LCL-filter grid-connected converter 


Table 1. The KP and Ki values 


Controller Kp Ki 
Active power controller 10 500 
Id controller 0.1 50 


Figure 6 depicts a comparing simulation results in Figure 6(a) the grid and injected active power, it 
was observed the times when the loads enter the work and how the grid is prepared for them, but at the time 
of the need to inject power. While, Figure 6(b) shows the injected active power and the performance of the 
controller was excellent and it gave the parameters: overshoot=1.08%, settling time 0.03 sec, rise 
time=0.00035 sec. Figure 7 comparing simulation results in Figure 7(a) it shows the grid's reactive power and 
the injected reactive power Figure 7(b), as there is no power injection in this work and its value is zero. 
Figure 8 shows the simulation results of the dq-axis current control with PI controller with trial-and-error 
method Figure 8(a) the d-axis control current and Figure 8(b) the q-axis current control. 
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(a)_ (b) 
Figure 6. System active power, (a) grid active power and (b) the injected active power 
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The current on the d-axis (ig) is what affects the injected effective power (P), and also the q-axis 
current (ig) for injected reactive power (Q). Therefore, controlling the current is controlling the power, and 
this form is identical. So controlling current is controlling power. Figure 9 shows the currents after using the 
PI controller with trial and error method and how the current increases during the period of increasing high 
loads and decreases during the period when the loads are out of work which Figure 9(a) the filter current and 
Figure 9(b) the grid current. 

By analyzing the harmonic of the converter current and the filter current, it was detected that the 
improvement of the voltage is obviously noticeable, as the voltage drop has become much lower, furthermore 
to eliminating high harmonics and keeping the necessary fundamental harmonic. THD was reduced from 
10.29% to 1.67% and this last rate is permissible according to IEEE 929-2000 and IEEE 519-1992 
specifications, as THD is less than 5% as presented in Figure 10, where Figure 10(a) presents the THD of the 
converter current, and Figure 10(b) depicts the LCL filter current. 
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Figure 7. System reactive power, (a) grid reactive power and (b) the injected reactive power 
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Figure 8. Simulation results of the dq-axis current control with PI controller with trial-and-error (a) D- axis 
required current and (b) Q-axis required current 
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Figure 9. Currents after using the PI controller with trial and error (a) LCL-filter current and (b) grid current 
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Figure 10. The THD, (a) Converter current and (b) LCL filter current 


6.2. PI-PSO controller with multi-objectives fitness function 

Then several intelligent technologies, such as PSO technology, have been presented to increase the 
capabilities of classical PI parameter tuning procedures in order to improve PI tuning. ITSE is the PSO 
objective function employed in this study. ITSE is expressed by the components integrator, product, clock, 
math function and work space blocks. To determine the best controller parameters, the PSO method employs 
sub-M-file programs. The following factors are used in the optimization method where C1, C2 equal to (1.3), 
W is (0.9), number of iterations is (100) and number of particles is (30). Table 2 shows the obtained 
parameters of the tow PI controllers. The tendency of particles moving within the searching field to the 
global optimal solution is depicted in Figure 11, which tracks the corresponding particles of the global 
optimal solution (the shown particles are part of the total particles). 


Table 2. Optimum controller parameters 


Controller Kp Ki Fitness 
Active power controller 7.45 735.64 63.22 
Id controller 0.3305 55.8527 


Figure 12 shows the controller performance, in Figure 12(a) the injected active power is depicted, 
Figure 12(b) shows the infected reactive power. Obviously, it performed admirably and provided the 
following parameters: overshoot=1.035%, settling time=0.015 sec, rise time=0.003 sec. This is an excellent 
result for a steady state error, but the overshoot is somewhat high. Figure 13 shows the simulation results of 
the dq-axis current control with PI-PSO controller with multi-objectives fitness function Figure 13(a) the d- 
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axis control current and Figure 13(b) the q-axis current control. Figure 14 shows a total harmonic distortion; 
Figure 14(a) shows the converter current THD, and Figure 14(b) shows the LCL filter current THD. These 
show that the THD was reduced from 10.23% to 1.57%. Finally, Figure 15 illustrates a performance 
comparison between the proposed controller and the traditional controller. 
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Figure 11. Particles tend to converge on the global optimum 
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Figure 12. Controller performance (a) injected active current and (b) the injected reactive power 
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Figure 13. Simulation results of the dq-axis current control with PI-PSO controller (a) D-axis control current 
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Figure 14. The injected current THD (a) converter current and (b) LCL filter current 
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Figure 15. Performance comparison of the injected active power between the proposed and the traditional 


controllers 
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7. CONCLUSION 

LCL filter is selected for usage in renewable energies applications, because of its high features. PI 
controller with trial-and-error method give 1.08%, 0.03 sec and 0.00035 sec for overshoot, settling time, rise 
time respectively and improving THD from 10.29% to 1.67% and PI-PSO controller with ITSE objective 
function give (63.2203996), (1.035%), ( 0.015 ) sec, ( 0.003 ) sec and from 10.23% to 1.57% for overshoot, 
settling time, rise time and improving THD respectively ,very exceptional in filtering from unwelcome 
harmonics and reducing the THD to less than 5%, according to the international standard for IEEE 929-2000 
and IEEE 519-1992, in addition to that it is light in weight and low in cost compared to other types. 
Depending on the results obtained by working in this paper the PI controller results are great especially in 
linear applications. Since the control of injected power tight and balanced, it does not change as demined 
even after operating the load in the grid. 
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